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We consider contributions to the rare decay K + — ► n + vT> which become non- 
. local at the charm scale. Compared to the leading term, such amplitudes 

are suppressed by powers of m 2 K /m%. and could potentially give corrections 
'"q! a ^ l eve l °f 15%. We compute the leading coefficients of the subleading 

dimension eight operators in the effective theory below the charm mass. The 
matrix elements of these operators cannot all be calculated from first princi- 
. pies and some must be modeled. We find that these contributions are likely 

J> , to be small, but the estimate is sufficiently uncertain that the result may be 

as large as the existing theoretical uncertainty from other sources. 
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The search for New Physics relies on experimentally accessible quantities whose Standard 
Model values can be predicted accurately and reliably. This task is often complicated by 
nonperturbative hadronic physics, especially when one is interested in the parameters of 
the Cabibbo-Kobayashi-Maskawa (CKM) matrix. To make progress, it is important to find 
processes where symmetry can be used to treat low energy QCD effects in a controlled and 
systematic way. One of these is the rare decay K + — > 7i + uu. This process is an example 
of a neutral current AS = 1 transition, which in the Standard Model can occur only via 
one-loop diagrams. 

The leading contributions to the effective Hamiltonian for this decay are given by 

w e// = % 9 a 2n E (y: s v u x{ Xt ) + v; s v cd x l NL ( Xc )) sY(l - 7 5 )^7,(1 - l>i , (i) 

\J I Z7T Sill K3W I 

where the index / = e,/z, r denotes the lepton flavor. The coefficient X(x t ) arises from 
the top quark loop and is independent of lepton flavor. It is dominated by calculable high 
energy physics, and has been computed to 0(a s ) |L|. Because it grows as mf, it is large 
and gives the leading contribution to the decay rate. If this were the sole contribution, 
the measurement of K + — > 7i + uu would yield a direct determination of the combination of 
CKM parameters |V^*V^| @. However, due to the smallness of V t * s V t d compared to V* s V c d, 
the charm contribution contained in the coefficient function X l NL (x c ) is significant as well. 
These terms have been calculated to next-to- leading logarithmic order ||. An important 
source of error in the calculation comes from the uncertainty in the charm quark mass, on 
which X l NL (x c ) depends. 

An important feature of the calculation is the fact that the hadronic matrix element 
(7r + |s7 I/ (l — rf)d\K + ) is related via isospin to the matrix element (7r°|s7 I/ (l — rf)d\K + ) re- 
sponsible for K + — > 7i°e + u. This largely eliminates the uncertainty due to nonperturbative 
QCD, up to small isospin breaking effects ||]. However, there remain long distance contri- 
butions associated with penguin diagrams containing up quarks which can lead to on-shell 
intermediate states. Some of these have been estimated in chiral perturbation theory and 
found to be small [0]. The perturbative contribution from virtual up quarks is tiny, since it 
is suppressed compared to the charm contribution by w? u jw? c . 

Summed over neutrino species, the branching fraction for K + — > n + uu is given by 



B{K + -> ti + vv) = k a 
with A = sin^c ~ 0.22 and 



'^AV,) + (^J2X^x c ) + ^X(x t ) 



(2) 



3a 2 B(K+ -»■ 7T°e + u) x8 
2tt z sin <d w 

Here & = V* s V ic i, and r K + absorbs isospin breaking corrections to the relationship between 
the decays K + — > 7r°e + z/ and K + — > tt + vv calculated in Ref. Q. In terms of the Wolfenstein 
parameterization of the CKM matrix |J, the branching ratio may be written as 

B(K + -> = 4.11 x 10- 11 • A 4 X 2 (x t )- \(af]) 2 + (p - pf] , (4) 
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FIG. 1. Penguin and box diagrams responsible for K + — > tt + vu. 



with 

a = (1- X 2 /2)~ 2 and p = l + 5 c , 

where 5 C absorbs the charm contribution. A measurement of the branching ratio then con- 
strains the parameters p and fj, which are equal to the Wolfenstein parameters p and rj up to 
known corrections of 0(X 2 ). The Alternate Gradient Synchrotron (AGS) experiment E949 
at Brookhaven and the CKM collaboration at Fermilab propose to obtain measurements 
of the branching ratio for K + — > 7i + uu at the level of 30% and 10%, respectively. The 
Brookhaven experiment is the successor to AGS-E787, which saw one event in this chan- 
nel [?] . These experimental prospects then fix the goal for the accuracy of the theoretical 
prediction at less than 10%. 

The leading source of theoretical uncertainty is associated with the charm contribution. 
Calculations at next-to-leading order in QCD yield S c = 0.40 ± 0.07, where the error is due 
primarily to the uncertainty in the charm mass 0. The errors from uncomputed terms of 
order a 2 (m c ) are expected to be small. However, the computation of the charm contribution 
relies on an operator product expansion which is simultaneously a series in a s and an expan- 
sion in higher dimension operators suppressed by powers of m c . The operators which are of 
higher order in the l/m c expansion reflect the fact that the penguin loop becomes nonlocal 
at the relatively low scale m c . One might expect the leading correction from higher order 
terms to give a contribution to 5 C of relative size m 2 K /m 2 c ~ 15%, large enough to affect in 
a noticeable way the extraction of p and fj from the decay rate. It is important either to 
verify or to exclude the presence of new terms of such a magnitude. 

In this note we will study the contributions of dimension eight operators to the decay 
K + — > 7i + uu. We estimate the correction to 5 C and comment on the uncertainty induced. 
After discussing the relevant power counting, we present the calculation of the operator 
coefficients and an estimation of the correction to the decay rate. We will find a small 
contribution, but one that need not be negligible. 

The decay K + — > n^vv proceeds via the loop processes shown in Fig. [I], which mediate 
the quark level transition s — *> dvv. These diagrams contain both short distance and long 
distance effects, which we separate by computing the effective Hamiltonian density 7i e s at a 
low scale p < 1 GeV. The effective Hamiltonian will receive corrections from the charm and 
top quarks, both of which have been integrated out of the theory, and from highly virtual 
up quarks. Soft up quarks remain in the theory, and are responsible for long distance 
corrections. 

We construct the effective Hamiltonian with an operator product expansion. At leading 
order, the operator in 7i e ff which contributes to the decay is of dimension six, 
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0(6) = - 7 5 )c/z/ 7 ^(l - 7 5 )i/. 



(5) 



The t quark contribution to the coefficient of this operator is obtained by evaluating the 
diagrams in Fig. [I] at the scale fi = My/ pa m t and matching on to the effective theory below 
this scale. At the same time, the W and Z are integrated out of the theory, producing four- 
fermion operators involving up and charm quarks as well. The charm contribution to the 
operator is then obtained by evaluating the diagrams contributing to the decay at /i = m c . 
These diagrams look like those in Fig. |Tj, but with the W and Z propagators replaced by 
local interactions. 

Dimensional analysis indicates that the coefficient of the dimension six operator O® 
scales as 1/Myy. The diagrams in Fig. [l| are quadratically divergent in the effective theory 
below Mw, and scale as A 2 /M v/ , where A ~ My/ is an ultraviolet cutoff. The Glashow- 
Iliopoulos-Maiani (GIM) mechanism ensures that this leading divergence cancels, since it is 
independent of the mass m q of the virtual quark. The consequence is that the coefficient 
of O® actually scales as rn^/M^. In terms of the Wolfenstein parameter A, the top coef- 
ficient has strength X 5 m 2 /M v/ and the charm coefficient has strength Xm^/M^. The top 
contribution is significant because of the large top mass, since X 4 m 2 /m 2 is of order 10. 

For the purpose of power counting, the operators of dimension eight scale as 

O® ~ M 2 K 0^ , (6) 

appearing with generic coefficient 0(8). Dimensionally, 0(8) is proportional to The 
top contribution to 0( 8 )0® is suppressed by M^/m 2 relative to its contribution to 0®, 
leading to an overall strength of order X 5 M] ( /M V/ . The corresponding suppression for charm 
is only M^/m 2 , so the overall contribution of charm to (7(8)0® scales as XM^/My/. Note 
that there is now no relative enhancement from the large top mass, so the top contribution 
to 0( 8 ) is suppressed relative to that of charm by A 4 and can be neglected. Furthermore, 
the contributions in question are independent of m q , so they cancel by the GIM mechanism 
when the up contribution is included. 

However, the GIM cancellation is manifest in Feynman diagrams only for contributions 
which are perturbatively calculable. The long distance contributions involving soft up quarks 
will differ by factors of order one from their perturbative representations. For these parts of 
the diagrams, which scale as 1/M^, the GIM cancellation is ineffective. Such long distance 
contributions have been considered elsewhere 0, and estimated to be small. The GIM can- 
cellation is also spoiled by logarithmic contributions proportional to (M^/M^) In (ml/M^). 
Such terms may be generated by the running of 7i e s between the scale m c and the low energy 
scale fj, < 1 GeV. This is not a large logarithm, numerically, but it allows us nonetheless to 
identify a GIM violating contribution to 7i e fr- This term, which is generated by intermediate 
up quarks as shown in Fig. ||], is of the same power-counting size as the long distance con- 
tribution. But because the perturbative description of the long-distance part is inaccurate, 
there is no reason to expect the GIM cancellation to be restored when it is included. 

The purpose of this paper is to compute the corrections to K + — > ix+vv of order 
(Mk/Mw) In (ml/M^). These contributions are well defined, and it is important, in light of 
the experimental situation discussed above, to determine whether they introduce a theoreti- 
cal uncertainty at a level competitive with the uncertainty due to m c . Note that pure power 
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FIG. 2. Diagrams leading to operators of dimension eight. 

counting arguments permit a relative contribution to 5 C of the order of (wf^/m 2 ) ln(m 2 /// 2 ), 
which could be as large as 20%, depending on the value chosen for the hadronic scale \x. 

We will study the effective Hamiltonian of dimension eight operators, at leading order 
in ol s . This Hamiltonian receives logarithmically enhanced contributions from the up quark 
loops in Fig. |^. We also must consider the matching corrections at the scale m c « m T , when 
the tau lepton is integrated out of the theory. Because the matching function F(m c /m T ) can- 
not be approximated by an expansion in m c /m T , the combination [F(m c /m T ) — F(m u /m T )} 
is a GIM violating finite matching correction which also must be included. 

The effective Hamiltonian density at the scale /i takes the form 

7ieff = E c ?0')Qi(A*), (7) 

l,i 

where / denotes lepton flavor. As it turns out, there will be two dimension eight operators 
generated in the theory below m c , 

0[ = sY(l - l 5 )d {idf [v llv {l - 7 >, , 

0\ = S7"(l - j 5 )(iD) 2 dv llv (l - 7 >, + 2s 7 "(l - j 5 )(iD»)du a »{l ~ 7 5 )(^>* 

+ s 7 *(l-7 5 )d^(l- 7 5 )(zd)V (8) 

The first of these operators does not receive any logarithmic QCD corrections below the 
scale m c , because it is proportional to a current which is partially conserved. The second 
does, but we will not include higher order corrections of relative order a s ln(m c //i). Note 
that this is not inconsistent with resumming terms of order a" ln n (M^/m c ). 

The operators 0[ 2 are generated by the diagrams in Fig. 0. In principle, one might have 
expected the diagrams in Fig. [3] to generate additional operators with a gluon field strength, 
such as 

0\ = sYcr ap G ap (l - 7 5 )rf iry„(l - 7 >z • (9) 

However, it turns out that contributions to all such operators cancel. 
The operators of dimension six that will induce 0[ 2 in TC e g are 

4 = s7"(l-7 5 )rfM 7! ,(l-7 5 )u, 

5 = ^(l- 7 5 ) M ^(l-7 5 )rf, 

0[ = u lv (-| sin 2 6 W (1 + 7 5 ) + (1 - | sin 2 6 W )(1 - 7 5 )) u v a v {l - 7 5 )^ , 

1 7 = - S1 V {\ -j 5 )u v llv {l 

1 8 = uY(l - l 5 )dhu(l ~ 7>z • (10) 



4 



V 



V 



V 



V 



FIG. 3. Diagrams which could lead to an operator with a gluon field strength. 

The operator 0\ comes from virtual Z exchange, the others from W exchange. The renor- 
malization group equations for 0[ and 0\ are 

^ = = 7 2 ^ + E^. (11) 

The anomalous dimensions 71, 7^ and 72 are of order one. The matrix 72 j is of order a s and 
comes from QCD running below m c ; it will not be included in our analysis. 

Computing the diagrams in Fig. |2| and solving the renormalization group equations, we 
find the coefficients at the scale /1, 



C ^' T = F772~(l ~ I sin2 ®w)G(a s ) log(/i/m c ) 



where 



„ ./,.- _ C I. . ; , , 2 , 

6M 2 , 

<^ = -^logWm c ), (12) 
d = -J£- } (ml/ml) , 



G F a 

V2 2tT Sill uyy 

r/ s f 6x -2 \, 4a; 

= 7 7v7 " 2 logx - -. — , (13) 



G(a s 



/ u (x-1) 2 

/ 1 \ \ -6/25 / / \ \ -6/23 / / v \ 12/25 / / x \ 12/23 

I a s (m c ) \ I ajjrib) \ I a s {m c ) \ I a s {m b ) \ 

\ a s (m b ) J \a s (M w ) J \a s (m b ) J \ a s (M w ) J 



and we have used V* s V u d ~ —V* s V c d- Taking the values m c = 1.3 GeV, mj = 4.5 GeV, 
A^jg = 0.35 GeV and // = m c /2, we find 

C^ r = 0.05 • c /M^ , = 0.69 • c /M^ , C T 2 = 0.28 • c /M^ . (14) 

By comparison, the coefficient of the leading charm contribution in Eq. (|l|) is given by 
X l NL (x c ) c , which is 4.0m 2 • cq/M^ for I = e,/z and 2.7m 2 • cq/M^ for Z = r. 
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To compute the contribution to the decay rate, we also need the matrix elements 
(tt + i>iVi \ 0\ 2 \K + ). The leading relative corrections come from the interference of 0\ 2 with 
and depend on the ratios 

Re MRSj Ktt+^^I |A-+)*(7r+^|Q(6)|i^+) 

/d[P.S.] |(7T+Z/Z7|0( 6 )|ir+)| 2 ' 1 J 

The matrix element of the operator 0\ is easy to calculate, since it depends only on the 
lepton momenta. The leptons are treated perturbatively, so the hadronic dependence of the 
matrix element of 0\ is the same as that of (6) . We then find 

Ri = (fa + Puf) = (180 MeV) 2 . (16) 

Unfortunately, the matrix element of 2 cannot be calculated analytically, since it in- 
volves the gluon field through the covariant derivative acting on the down quark. We are 
forced to rely instead on model dependent estimates, which are notoriously unreliable. One 
ansatz would be to take 

(n + uu\0 2 \K + ), « p 2 (it + uu\0^\K + )^ (17) 

or R 2 ~ fi 2 ~ (650 MeV) 2 . Another would be to neglect the gluon field and model the 
matrix element as 

(u + uu\0 2 \K + ) = fa+p P ) 2 (-k + up\Q^\K+), (18) 

in which case 

R2 ~ {fa + Pd) 2 ) = (340 MeV) 2 . (19) 

Of course, neither of these guesses need be correct within better than an order of magnitude. 
Fortunately, lattice QCD methods are advancing quickly, to the point that a true unquenched 
lattice calculation of this matrix element may soon be feasible. For now, we will take these 
two crude guesses to bracket roughly the actual value of R 2 . 

We now write the branching fraction for K + — > 7i + uu as in Eq. (^), with 

po = 1 + 4(1 + h) , (20) 

where 5§ is the new term which we are computing. Summed over lepton species, the contri- 
bution of charm at dimension six is given by 

6c = A°X(xt) = 3A* £*^(*e) " A2x{xt) ■ (21) 

A next to leading order analysis yield Po — 0.42 ± 0.06, where the error arises in large part 
from the uncertainty in the charm quark mass 0. This value of Po gives S c = 0.40 ± 0.07, 
where we use X(x t ) = 1.53 ± 0.01 and A = 0.83 ± 0.06. The fractional correction due to 
dimension eight operators is then 
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s * = ^ E ( c i R i + = 4 1} + <g° • (22) 

The first term, for which the matrix element is calculable, is negligible in size: with our 
choice of inputs, 5^ = 5.6 x 10 -5 . The second, highly uncertain, term is much bigger, with 
5^ between 1% and 5% for our adopted range for R 2 . On the one hand, even 5$ as large as 
5% is somewhat below the existing uncertainty on 5 C from the value of m c . On the other, if 
our "upper limit" on the matrix element of O l 2 were too small by even a factor of two, which 
need not be unlikely, these contributions would have a significant effect on the extraction of 
CKM parameters from the branching fraction. 

We have made a number of approximations in obtaining these results. A potentially 
important one, within the perturbative calculation, is that we have neglected QCD running 
below m c . We could include these QCD corrections for 0% simply by incorporating the 
known running of the coefficients C 4 and C 5 ; doing so decreases C\ by a factor of two. 
However, the running of 2 is not equally trivial, since 2 itself is renormalized in QCD. In 
view of the large uncertainty in the matrix element of 2 , including these QCD corrections 
would not at this time increase the reliability of our prediction. 

Of course, the key uncertainty arises not from QCD perturbation theory but from the 
actual value of (tt + uu\ 2 \K + ). Only a realistic lattice computation will settle the matter. 
We would argue, in fact, that such a calculation is really required for one to be confident 
that the effects we have considered do not spoil the extraction of CKM matrix elements 
from the proposed experiments on K + — > 7i + uu. This is not the only case where higher 
dimensional operators can play an interesting role in kaon decays 

In summary, we have computed the dominant contribution to the coefficients of dimension 
eight operators contributing to the decay K + — > tt + uu. Our best estimate is that this 
represents a correction of no more than 5% to the leading charm contribution to the decay. 
However, our ignorance of relevant hadronic matrix elements leaves open the possibility 
that these contributions could represent an uncertainty as large as or larger than that due 
to the charm quark mass. A lattice calculation of nonperturbative corrections, and to a 
lesser extent the inclusion of perturbative QCD corrections below the charm scale, will be 
indispensable to reducing this uncertainty before the planned experiments begin to take 
their data. 
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